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■ It is shown that the A^-dark soliton solutions of the integrable discrete nonlinear Schrodinger 

(IDNLS) equation are given in terms of the Casorati determinant. The conditions for reduction, com- 
plex conjugacy and regularity for the Casorati determinant solution are also given explicitly. The 
relationship between the IDNLS and the relativistic Toda lattice is discussed. 
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1. Introduction 

The nonlinear Schrodinger (NLS) equation, 



(-sq ; tyt = Yxc + o#|>, (1-1) 

in 

is one of most important soliton equations in mathematics and physics. The study of discrete analogues 
of the NLS equation has received considerable attention recently from both physical and mathematical 



point of view. 1 - 2 The integrable discrete nonlinear Schrodinger (IDNLS) equation is given by 

d\\f n ,o, 

i— - =\|f„+i +V»-l + a¥n (Vn+l +Vn-l)- (1-2) 
at 

The IDNLS equation was originally derived by Ablowitz and Ladik using the Lax formulation. 1 ' 3 ' 4 
Thus the IDNLS equation is often called the Ablowitz-Ladik lattice. 



X 

The IDNLS equation can be bilinearized into 



D t G n F n — G„ + iF„_i —F n+ \G n -\ , (1.3) 

FZ + aG n G* n = F n+{ F n _y, (1.4) 

through the dependent variable transformation \|/„ = (i)"G n /F n where F n is a real function, G n is a 
complex function and G* is its complex conjugate. The Hirota D-operator D t is defined by 

D?f-g = {d t -Z t ,) n f{t)g{t% =t . 

With this bilinear forms the IDNLS equation in the case of a = 1 can be solved by the Hirota bilinear 
method, i.e. Af-bright soliton solutions are obtained. 5 "" 10 
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It is well known that the NLS equation with defocusing parameter (a < 0) has dark soliton solu- 
tions. It is easily confirmed that there is a stationary dark soliton solution 



by the direct substitution of an ansatz of stationary solution into eq.(1.2). Understanding various prop- 
erties of dark solitons in discrete lattices is important from physical point of view. 11-13 There are 
analytical results of the IDNLS equation under non-vanishing boundary condition using the perturba- 
tive Hirota method by Narita 14 and inverse scattering theory by Vekslerchik and Konotop. 15 However, 
the detailed anayisis of the solution and understanding the role of a lattice space parameter in the N- 
dark soliton solutions are missing in their results. In ref.15, a determinant expression of dark soliton 
solutions is given, however for their boundary condition there is no carrier wave of the dark soliton. 
Besides this is a Grammian type determinant solution. The Casorati determinant solution which is 
well-known in discrete integrable systems is missing. 

The relationship between the IDNLS and relativistic Toda lattice (RTL) equation was also pointed 
out by several authors. 101617 An interesting open problem was given in conclusion in ref.10: The 
2-component Casorati determinant solution of the IDNLS equation and RTL equation in ref.10 is 
different from the single-component Casorati determinant solution of the RTL equation which was 
derived in ref.18. Is there any relationship between the 2-component Casorati determinant solution 
and single-component Casorati determinant solution? 

Our goal in this paper is to construct Casorati determinant form of Af-dark soliton solutions with a 
lattice space parameter, analyze the detail of behaviour of Af-dark soliton solutions and give an answer 
of the above open problem. In this article, we give an explicit formula of the dark soliton solutions. It is 
known that the dark soliton solutions for NLS are written in terms of Wronski determinants. We show 
that the solution for IDNLS is given in the Casorati determinant form which is a discrete analogue of 
the Wronskian. 

This paper is organized as follows. In the next section, we discuss about the gauge transformation 
and bilinearization of IDNLS equation. In §3, we give bilinear identities for Casorati determinant both 
for the Backlund transformation of Toda lattice and the discrete 2-dimensional Toda lattice equation. 
In §4, we perform a reduction of the Casorati determinant which derives the bilinear form for IDNLS 
given in §2. The complex conjugate condition is considered in §5. In §6, we give examples of dark 
soliton solutions. We also clarify the relationship between the IDNLS and the relativistic Toda lattice 
in the final section. 

2. Gauge Transformation of IDNLS 

First of all, we should examine the transformation of IDNLS equation by gauge. Applying the 
gauge transformation, 



\|/„ = exp(-2ifsech 2 (fc)) 



tanh(fc) 



tanh(kn) 



(1.5) 




\|/„ = w„Aexp(i(9ra — CO?)), 
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eq. (1.2) is transformed as 

.dun _ _ (0u ^ _|_ -\-aA 2 \u n \ 2 )(u n+ i exp(i6) + exp(— i6)) , 
dt 

thus we can scale | oc| by using the amplitude A. To obtain the dark soliton solutions, we set a = 
— (a 2 — l)/a 2 , (a > 1). We also scale t by t — > a 2 t, then the IDNLS eq. (1.2) can be bilinearized into 



(ID, + co)G n • F„ = G„+iF„_! exp(i0) + G„_iF„ +1 exp(-i6) , 



(2.1) 
(2.2) 



through the dependent variable transformation, 

G n 



V B = — exp(i(6n-aw)), 

where F„ is a real function, G„ is a complex function and G* is its complex conjugate. 

The above gauge factor exp(i(6« — aw)) is nothing but the carrier wave of the dark soliton, and 
essentially u n = G n /F n gives the envelope of the soliton solution. In the context of NLS, in order to 
construct physical solutions, it is more convenient to consider the above bilinear form than the one in 
eqs. (1.3) and (1.4). Using the bilinear forms, eqs. (2.1) and (2.2), we present the Casorati determinant 
form of dark soliton solutions of the IDNLS equation. Instead of applying the Hirota bilinear method 
directly to eqs. (2.1) and (2.2), we rather start from the algebraic identities for Casorati determinant 
and derive the bilinear form of IDNLS equation by using the reduction technique. 

3. Bilinear Identities for Casorati Determinant 

3.1 Backlund transformation ofToda lattice 

We have the bilinear forms for Backlund Transformation (BT) of Toda lattice (TL) equation, 



(aD x - 1)t b+ i (k+ 1,1) -% n (k,l) +% n (k+ 1, l)x n+l (k,l) 

(bD y - l)T B _i (*, I + 1) ■ T„(fc, /) + X n (k, I + \)% n -l(k, I) 

where a and b are constants. 

Consider the following Casorati determinant solution, 



0, 
0, 



(3.1) 
(3.2) 



T«(M) = 



q>W(M) <p! n+1) (M) - cpf^ 
<#>(*,/) <P? +1) (M) ••• cp^- 1 ' 



(*,/) 
(*,/) 

(*,/) 



(3.3) 



where cp- s are arbitrary functions of two continuous independent variables, x and y, and two discrete 
ones, k and /, which satisfy the dispersion relations, 

d x ^\k,l) = ^ +l \k,l), (3.4) 
d y <vt\k,l) = <Q { r l \k,l), (3.5) 



(»+i) 



A / c P( w (^/)=( P r i; (^o 



(n-l). 



(3.6) 
(3.7) 
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Here A k and A/ are the backward difference operators with difference intervals a and b, defined by 

/(*,/)-/(*-!,/) 



Akf(k,l) = 



Aif(k,l) 



f(k,l)-f(k,l-l) 



(3.8) 
(3.9) 



For simplicity we introduce a convenient notation, 

q>S Bl) (*i,/i) <? { " 2 \k 2 ,h) 
9?°(*i,/i) <? { 2 2 \k 2 ,l 2 ) 



>2) 



("at)/ 



(3.10) 



In this notation, the solution for BT of TL, T n (k,l) in eq. (3.3), is rewritten as 

T„(M) = K/'" + 1 M'"' ,« + #-!*,/!> (3-11) 
or suppressing the index k and I, we simply write as 

l n (k,l) = \n,n + 1, ■ ■ • ,n+N— 1|. 

We show that the above x n (k,l) actually satisfies the bilinear eqs. (3.1) and (3.2) by using the 
Laplace expansion technique. 19 20 At first we investigate the difference formula for x. From eq. (3.6) 
we have 



a(p! w+i) (k +i,o= q>r (* + 1 , o - q>r (*> o • 

Noticing this relation, we get 

T„(&+1,/) = |n£ + i,/,n+ l fe+ i )/ ,n + 2 fe+ i ) /,--- ,n+N- 

= \nk,l,n+lk+i,l,n + 2 k+l j,- ■■ ,n + N-l k +i,i\ 

where we have subtracted the 2nd column multiplied by a from the 1st column, 

= \n u ,n+\ u ,n + 2 k+u ,--- ,n+N— 

where we have subtracted the 3rd column multiplied by a from the 2nd column, 



(»), 



that is, 



= K,/,n + !*;,/,■•• ,n + N-2 k j,n + N- l k+1 



T n (k+l,l) = |nfe,/, « + !*,/'■"■ ,n + N-2 k j,n+N- 



(3.12) 



(3.13) 



Moreover in eq. (3.13), multiplying the iV-th column by a and adding the (AT — l)-th column to the 
N-th column, we get 



az n (k+l,l) = \n k j,n+\ k j,--- ,n+N -2 k j,n+N -2 k+u \ 



(3.14) 
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Differentiating eq.(3.14) withx and using eq.(3.4) we obtain 

ad x % n (k+l,l) = |n*,/, n + lM'"" ,n+N -3 k j,n + N - l kjh n + N- 2 k+1 j\ 
+ \n k j,n+l u ,--- ,n+N-3 k j,n+N-2 k j,n+N-l k+u \ 
= \n kh n + \ kh --- ,n+N-\ u n + N-\ kh n + N-2 k+ ij\+x n {k+\,l). (3.15) 

We have also 

d x x n {k,l) = \n k j,n+l k j,--- ,n+N -2 k j,n+N k j\. (3.16) 

In short, we write 

% n (k+l,l) = \n,n + l,--- ,n+N-2,n+N-l k+1 \, (3.17) 

a% n (k+\,l) = \n,n + l,--- ,n + N -2,n + N -2 k+i \ , (3.18) 

(ad x -l)i n (k+l,l) = \n,n + l,--- ,n + N-3,n + N- \,n + N-2 k+l \ , (3.19) 

d x T„(k,l) = \n,n+l,--- ,n + N-2,n + N\. (3.20) 
Let us introduce an identity for 2N x 2N determinant, 
n+\ ■■■ n+N-2 \ n+N \ n+N-\ k+x i i n+N-l 

i n+N i «+A^-lyt + i i n n+l ■■■ n+N-2 i n+N-l 

Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear identity for 
determinants, 

|n + 1,- - ■ ,n + N-2,n + N,n + N-l k+ i \ x \n,n + 1,- • ■ ,n + N -2,n + N - l\ 
\n+l,--- ,n+N-2,n+N-l,n+N-l k+l \ X \n,n+l,--- ,n+N-2,n+N\ 
- |n+ 1, ■ - ■ ,n + N -2,n + N,n + N-l\ x \n,n+ 1, ■ ■ ■ ,n+N-2,n + N- l k+ i \ =0, (3.21) 

which is rewritten by using eqs. (3.17)-(3.20), into the differential bilinear equation, 

(ad x - 1)t b+ i(*+ 1,/) x x n (k,l)-ax n+l (k+ l,l)d x x n (k,l) +i n+1 (k,l)l n (k+ 1,1) = 0. 

This equation is equal to eq. (3.1). 
Similarly we obtain 

(bd y - l)x„_i (k, l + l)xT n (k,l)- bi n -i(k, I + l)dyi n (k, I) + z n -i(k, l)x n (k,l + l) = 0, 

because the role of y, I and b is parallel to that of x, k and a except that the ordering of index n is 
reversed. The above equation is equal to eq. (3.2). Hence we have proved that the Casorati determinant 
T n (k, I) in eq. (3.3) gives the solution for the BT of TL equation. 



0. 
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3.2 Discrete 2 -dimensional Toda lattice 

Using the Casoratian technique, we show that the bilinear form of discrete 2-dimensional Toda 
lattice (D2DTL) equation, 

x n {k+\,l + \)x n {k,l)-x n {k+\,l)x n (k,l + \) 

= ab{i n (k + 1 , / + 1 )t„ (*, / ) - t b+ i (* + 1 , / )t„_ i (*, / + 1 ) ) , (3.22) 

is also satisfied by the Casorati determinant. 

Let us further examine the difference formula for the Casorati determinant x„ in eq. (3.3). Similarly 
to eqs. (3.17) and (3.18), we have 

x n (k,l + l) = \n i+1 ,n + l,--- ,n + N-2,n + N-l\, (3.23) 

b%„(k,l + l) = |n + l/+i,n+l,-- - ,n + N -2,n + N - \\. (3.24) 

Thus T n (k + 1,1 + 1) is given by 

% n (k+l,l + l) = \n k+ i t i + i,n + \ k+ i,--- ,n + N -2 k +\,n+N -\ k +\\. (3.25) 

In the following we show that the shifts of index k is condensed into only the right-most column of the 
determinant. From eqs. (3.6) and (3.7), q>^ satisfies 

AAcpJ n) =cpJ' l) , (3.26) 

that is, 

q>J B) (*,/)- <P, W (*-l,0- q>l B) (*,/-!) + q>l B) (* - 1 , / - 1 ) = ^cp| n) (*,/), (3-27) 
which is rewritten as 

(l-^)(p, W (^+l,/ + l) = (pW^^+ij + tpW^+i,/)-^^,/) 

= (pJ' l) (A:,Z + l)+a(pJ' l+1) (^+l,/). (3.28) 

Therefore multiplying the both hand sides of eq. (3.25) by (1 — ab) and rewriting the 1st column of 
the determinant by the use of eq. (3.28), we obtain 

(l-ab)T n (k+\,l + \) 
= \n k j + i,n+\ k+i ,--- ,n+N-2 k+u n+N-\ k+ i\ 

+ a\n + \ k+u ,n + \ k+u --- ,n + N-2 k+i ,n + N-\ k+ i\. (3.29) 
The second term in r.h.s. vanishes because its 1st and 2nd columns are the same. So we get 

(\-ab)l n (k+\,l + \) = \ni +u n + \ k+1 ,--- ,n + N-2 k+l ,n+N-\ k+l \. (3.30) 

In the above determinant, by subtracting the (j + l)-th column multiplied by a from the j'-th column 

for j = 2,3, ■ ■ ■ ,A^— 1, z n (k+ 1,1 + 1) is given by 

(\-ab)% n (k+\,l + \) = |n m , «+!,■■■ ,n + N -2,n + N - l k+1 \. (3.31) 
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As is shown, even if the two variables k and I are shifted, x n is also expressed in the determinant 
form whose columns are almost unchanged and only edges are varied. Thus we can use the Laplace 
expansion technique. 

We consider an identity for 2N x 2N determinant, 

m+i i n + 1 ■■■ n + N-2 n + N-\ k +i \ n\ i n + N-l 

n, +1 : in: n + l ••• n + N-2 i n + N-l 

By the Laplace expansion, we have 

|n/_i_i,n+ 1, - ■ ■ ,n + N — 2,n + N— ljt+i| x |n,n+ 1, ■ ■ ■ ,n+N — 2,n+N— 1| 
|n,n + l,--- ,n + N -2,n + N -l k+ i \ x |n/ +1 ,n + 1, • • • ,n+N — 2,n+N— 1| 
+ |n + l,--- ,n+/V-2,n+/V- l,n+/V- l fc+1 | x |n /+ i ,n,n + 1, • • • ,n + N-2\ =0.(3.32) 

Using eqs. (3.17), (3.18), (3.23), (3.24) and (3.31), we obtain from eq. (3.32), 

(1 — ab)% n (k + 1,1 + l)l n (k,l) —T n (k+ \,l)l n (k,l+ 1) 
+ai n+l {k+\,l)bx n -i{k,l + \) = 0, (3.33) 

which recovers eq. (3.22). This completes the proof that the Casorati determinant is indeed the solution 
for D2DTL equation. 

(n) 

For example, we can take q>- as 

4 n \W)=tf(l-a Pi )- k (l-bj^ exp(? 3i ) + <fi(\-aq i )- k (\-b^J exp^), (3.34) 

1 K 1 

t,i = piX + — y + ^o, X\i = q t x + — y + X\ iQ , 
Pi qt 

where qi and ^,o, r(,o are arbitrary constants which correspond to the wave numbers and phase 
parameters of solitons, respectively. It can easily be seen that cp-"' in eq. (3.34) actually satisfies the 
dispersion relations, eqs. (3.4)-(3.7). 

4. Reduction to IDNLS 

In this section, we give the reduction technique in order to derive the IDNLS equation from the 
system of BT of TL and D2DTL. In eq. (3.34), we apply a condition for the wave numbers, 

qi = -Pr, -• (4-1) 

1 - apt 

Then we get 

l-b l 

Pi = -q, 

and thus 



Pi = -qn q ~, (4.2) 

1 - aqi 



*=9h— (4-3) 

1 - api 1 - aqi 
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On this condition, q>- in eq. (3.34) satisfies 



^\ k+ l,l-l)= p f ^\k,l) 

1 G-Pi 



Hence we obtain 



T„ +2 (ft+l,/-l) = [Y[fi-—1L )x„(k,l). 
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(4.4) 



(4.5) 



By using eq. (4.5), the bilinear forms for BT of TL and D2DTL, eqs. (3.1), (3.2) and (3.22), are 
rewritten as 

(aD x - 1)t„+i (ft + 1, /) ■ T n (ft, l)+tn(k+h OVi (*, 0=0, (4.6) 

(bD y - l)T„+i (ft + 1, /) ■ T n (ft, I) + T n+2 (k + 1 , l)X n -l (*, 0=0, (4.7) 

x n+ i (ft + 1 , Z)x„_i (ft - 1 , 1) - x„+i (ft, Z)t„_i (ft, /) 

= aft(T n+ i (ft + 1 , Z)t b _i (ft - 1 , /) - T B (ft, l)z n (k, 0). (4-8) 

Here we may drop /-dependence, thus for simplicity we take I = hereafter. Equation (4.8) is nothing 
but the bilinear form of the discrete one-dimensional TL equation. 
Let us substitute 

x = iact, y = ibdt, (4.9) 
where c and d are constants. Then we have 

- id, = acd x + bdd y , (4.10) 
and from the bilinear eqs. (4.6) and (4.7), we get 

(-iA -c-d)% n+ i(k+ 1) -T„(ft) + cx„(ft+ l)T„+i(ft) + d% n+2 {k+ l)T„_i(ft) =0. (4.11) 
The Casorati determinant solution of eqs. (4.11) and (4.8) (with / = 0) is given by 

cpW(ft) cp? I+1 >(ft) ••• cpf 
<#>(*) cp? +1 '(ft) 



!„(*) 



.(»+!)/ 



(pW(ft) = tf(l -a Pi )-*exp&) + rf(l -a ?i )-"exp(Tl ; ), 

- . / bd\ „ ■( bd\ 

t,i = i I acp; + — 1 1 + ^o, r|i = i I + — 1 1 + r\ i0 . 



Here and are related by eq. (4.1). 
By defining 



/n = T„(0), g„=T„+i(l), A„=T„_i(-l), 



(4.12) 
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from eqs. (4.8) and (4.11), we obtain the bilinear form of IDNLS equation, 

(-LD t -c- d)g n -f n + dg n+ if n -i + cg n -if n+ i = 0, (4.13) 

(\D t -c- d)h n ■ f n + ch n+ if n -\ + dh n -if n+ i = , (4. 14) 

fn+lfn-l -fnfn = (ab - 1 )(/„/„- gnK) ■ (4-15) 

By using the dependent variable transformation, 

_ gn _ K 

u n — , V n — , 

Jn Jn 

the above bilinear equations are transformed into 

^dun = _^ C J rd ^ Un + ^ ab _^ ab _ i) MfjVn )(j Mn+1 +cw„_i), (4.16) 
dt 

_.dvn _ _^ c + d ^ Vn + ^ ab _ (ab-\)u n v n )(cv n+ {+dv n -{) . (4.17) 
dt 

5. Complex Conjugate Condition 

In order to take v„ to be the complex conjugate of u n , we need to restrict /„ as real and h n as the 
complex conjugate of g n up to gauge freedom. Moreover to take u n to be regular (i.e., u n not to diverge 
for real t), we need /„^0 for real t. In this section, we give the conditions for the complex conjugacy 
and regularity for the Casorati determinant solution. 

Firstly we take 

b = a , 

for simplicity, and next we take 

a > 1 . 

Then the conditions for complex conjugate are given as follows: 

Pi = a+ \/a 2 -\r u |r,-| = 1 , (5.1) 

exp(Tl ;0 ) = I ft ELl3h ) exp(-^o) , (5-2) 
yt=t <li -qk J 

d = c*, (5.3) 

where r, is a complex parameter of absolute value 1, and * means the complex conjugate. On the 
condition (5.1), eq. (4.1) turns to be 

1 

and eqs. (5.2) and (5.3) gives 

ex P (T l! )=(n^ Z f)ex P (-^), 

for real t. 
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Now we summarize the final result of the Casorati determinant solution for IDNLS equation. The 
bilinear form of IDNLS equation is given as 



(-LD, - C - C*)g n ■ f n + C*g n +lf n -l + Cg n -lf n +l = , 
(\D, -c- C*)h n ■ f n + ch n+ \f n -\ + C*/l„_i/„ + i = , 
fn+lfn-l — fnfn = (« 2 ~ l)(/n/« — gnK) , 



(5.4) 
(5.5) 
(5-6) 



where a > 1 and t is real. The Casorati determinant solution for the above bilinear equations is given 
as 



fn 



9r ; (o) 

#(0) cp? +1) (0) 



^(n+iV-l) 



cpW(o) q>r i; (o) • 
cp ( ; +1) (i) cp ( ; +2) (i) 

cp? +1) (l) cp? +2) (l) 



(»+!)/ 



cpj; +1) (i) <pjr 2) (i) 



(n+N-l) 
(n+N) 

q>i 



(0) 
(0) 

(0) 

(i) 



(5.7) 



qT^C-l) #(-1) 



CP? +A?) (1) 



cpr w) (i; 



cp ( i" + ^ 2) (-l) 
cp? +A, - 2) (-D 



(5.8) 



cpr i} (-i) #(-D - cpr-^-D 



(5.9) 



where 



ft 



ft r / 



= ia ( c Pi + — ) t + fyo, pi = a + Va 2 -lt 



1, 



exp(-£), 



where r,- is a complex parameter of absolute value 1 and is an arbitrary complex parameter. 
This Casorati determinant solution gives the A^-dark soliton solution for IDNLS equation, r,- and 
parametrize the wave number and phase constant of i-th soliton, respectively. 

After a straightforward and tedious calculation of expanding the determinants, /„, g n and h n can 
be expressed in an explicit way as follows: 



fn — F n g , g n — G n g 

where Q is the gauge factor defined by 



nr 2 ---r N 



N 



h n =H n g 



r\r 2 ---r N 



(5.10) 



N>i>j>l i=l 
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and F n , G„ and H n are given as 



N 



^=1 I 

M=0 \<h<i2<-<i M <N 



n- n* — 1 



M 



riKi 2 "e^ 



v=l 



N 

M=0 l<(' 1 <i 2 < - <iM<A' 
iV 

M=0 l</i<i2< - <iM<A' 



Pip - Piy 



n ... ,, _ i 

l</j<V<M P'tiPh 

n 



Pip-Piv 



n 

v=l 
M 

n 
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(5.11) 

\Pi v \ 2n ^ iv+ ^, (5-12) 
Kl^. (5-13) 



1<H<V<M P'fPi v ~ 

A proof of the above expression is given in appendix B B. Now it is clear that F n is real and H n = G*, 
thus v„ is the complex conjugate of u n up to multiplication factor. Moreover since F„ > 1, the regularity 
of u n is also satisfied. 

From the bilinear eqs. (5.4)-(5.6), and the gauge transformation (5.10), the above F n and G n satisfy 
the bilinear equations, 



Let us take 



(-iD, - c - c*)G n ■ F n + c*G„+iF„_i + cG n -iF n+l = , 
-F n F n = (a 2 - \){F n F n - G n G* n ) . 

c = exp(— i8), 



(5.14) 
(5.15) 



where is a real parameter. By the variable transformation, 



\|/„ = — ^exp(i(6n — (at)), 



that is, 



where 



(Va 2 - \) N 

\if n = u n exp(i(6« — (Ot)) , 

nr 2 ---r N 



(0 = exp(iS) +exp(-i0) =2cos0, 



we finally obtain the IDNLS equation, 

i dMf n 



a 2 dt 



X|/„+l+\|/„-l 



a 1 -! 



|Vb| (Vb+1 +Vn-l)- 



(5.16) 



Thus we have completed the proof that the A^-dark soliton solutions for the IDNLS equation are given 
in terms of the Casorati determinant. Here we remark that a parameter a is related to a lattice space 
parameter 8 = I /a (< 1). Using a lattice parameter 8, we can rewrite the IDNLS equation (5.16) into 



. dM[ n _ \|/„+i +\|/„_1 1 - 8 2 2 



|¥«l(¥n+l+¥n-l)- 



(5.17) 



dt 8 2 8 2 

We see that the IDNLS equation is defocusing whenever a lattice parameter 8 is smaller than 1. When 
8 = 1, the IDNLS equation reduces to a linear equation. When 8 > 1, i.e. a < 1, the IDNLS equation 
is focusing. In this case, above single component Casorati determinant solution does not satisfy the 
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Fig. 1. Plots of 1-dark soliton. The vertical axis is |\|/„|. ((a) a = 2,r = exp(i7t),9 = K, (b) a = 1.1, r = 
exp(i7i),8 = 71.) 



conditions for complex conjugacy (unless we admit singular solutions). However, our Casorati deter- 
minant form can be extended into the case of the focusing IDNLS equation. This is corresponding to 
the homoclinic orbit solution of the IDNLS equation. The detail will be discussed in elsewhere. 

6. Soliton Solutions 

In this section we give examples of dark soliton solutions for the IDNLS equation. By taking 
N=\in eqs. (5.11) and (5.12), we get the 1-soliton solution for the IDNLS eq. (5.16), 

\|/„ = -^exp(i(8ra-2fcos8)), 
F„ 

where is a real parameter and 

F n = l + Ipil^expfo , G n = \- J^lp^exp^i +&) , 

where 

e ■ f Pi exp(i8)\ e r~i — 7 i 

si =la — t^vH m+sio, Pi= « + V a -in, \n =1- 

Vexp(i6) pi / 

By rewriting as pi = /?exp(i0) and ri = rexp(i0), we get a slightly simpler expression, 
F n = 1 + H^exp^ + £*) , G n = 1 - -^|p|2»exp(^ + ^) , 

where 

^ = ia (/? + —) ? + ^o, = aexp(— i6) + ^/ a 2 — lr, = 1 . 

V P/ 

Here arg(r) parametrizes the wave number of soliton, and SR(^o) gives the phase constant which 
parametrizes the position of soliton. Figure 1 shows a stationary 1-dark (black) soliton. When a(> 1) 
is getting closer to 1, the width of a dip is getting wider. Figures 2 and 3 show travelling 1-dark (gray 
and black) solitons. 

Here we should comment about the difference of travelling velocity between black and gray soli- 
tons. In the continuous case, the NLS equation (1.1) is invariant under the gauge and Galilean trans- 
formations, 

t) — > \\f(x - 2kt,t) exp(i(fce - k 2 t)), 
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Fig. 2. Plots of travelling 1 -gray soliton. The vertical axis is |\|/„|. ((a) a = 2,r = exp(3i7t/4),8 = 57t/3.) 
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Fig. 3. Plots of travelling 1 -black soliton. The vertical axis is |\|/„|. (a = 2,r = exp(— 5i7t/3),6 = 5%/3.) 



thus the travelling velocity of any solution can be shifted by 2k by using the above transformation. 
If we normalize the freedom of Galilean transformation by requiring that the carrier wave disappears 
(i.e. k = 0) and \|/ represents only the envelope, then the velocity of the envelope soliton is if and 
only if the soliton is black. In the case of IDNLS, the situation is almost same. When = (i.e. the 
case of carrier-wave-less) or = 71 (i.e. \|/„ = (—1)" x (envelope soliton)), then the velocity of black 
soliton is and the travelling dark solitons are always gray solitons. We see this fact from Figs. 1 and 
3. 
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The 2-soliton solution is obtained by taking N = 2 in eqs. (5.11) and (5. 12). Similarly by rewriting 
as pi — > p,exp(i9) and r,- — ► r,exp(i6), we get an explicit expression of the 2-soliton solution for 
eq. (5.16), 

\|/„ = — exp(i(8«-2fcos8)), 



+ 



G n 



1 + l/nrexp^ + + l^l^expfe + %) 
p] ~ r:: |pip 2 | 2 "expfti+«+§2 + §5), 



= i- 

+ 



i 

pi 



Pi-Pi 



P\Pi 



1 



p\n Pi r i 



where 



£,i = ia ^pi + ~ S j t + ^0, Pi = aexp(-ie) + \J a 1 - \r u 



1. 



Here the carrier wave of the dark soliton solution is given by the exponential factor, exp(i(6« — 
It cos 6)). Figure 4 shows 2-dark soliton interaction. 

7. Concluding Remarks 

We have shown that the Af-dark soliton solutions of the IDNLS equation are given by the Casorati 
determinant. From the above derivation of the jty-dark soliton solutions for the IDNLS equation, we 
notice that there is correspondence between the Casorati determinant solutions of the IDNLS equation 
and the RTL equation. 18 The RTL equation is decomposed into three bilinear equations, i.e. the Toda 
lattice (TL) equation, the Backhand transformation (BT) for TL and the DTL equation (See appendix 
A A). The IDNLS equation is also decomposed into three bilinear equations, i.e. two types of the BT 
for TL and the DTL equation. 



This fact reminds us several works on the relationship between IDNLS and RTL 



10, 16, 17 



Moreover, 



our result is an answer to an open problem of ref.10, i.e., the Af-bright soliton solutions of IDNLS are 
written by the 2-component Casorati determinant in ref.10, and the A-dark soliton solution are written 
by the single-component Casorati determinant which is given in this paper and corresponds to the 
solution of RTL in ref.18. 

It is interesting to investigate A-dark soliton solutions of the discrete time IDNLS equation. It will 
be reported in the forthcoming paper. 
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Fig. 4. Plots of 2-dark soliton interaction. The vertical axis is |\|/„|. (a = 2,r\ = exp(3i7t/4),r2 = 
exp(i7t/4),6 = 5jt/4) 



Appendix A: Relativistic Toda Lattice 

In this appendix we briefly explain x-functions of the RTL equation. The RTL equation 

d 2 q„ 



dt 2 



, 1 u, , 1 . v exp( g „_i- g „) 
c c 1 + 4 exp(^r„_i 



- (l + -?n)(l + -?i.+i)- 1 ; r, (A-l) 

c c l + -^exp(q n -q n+ i) 

where g„ is the coordinates of «-th lattice point and c is the light speed, was introduced and studied by 

Ruijsenaars. 21 The RTL equation (A-l) is transformed into the three bilinear equations, 

D 2 x f n -f„ = 2(g ngn -f 2 ), 

(aD x - \)f n -fn-\ + gn-lgn = 0, (A-2) 
gn-lgn+l ~ f„= a 2 (fn+lfn-l ~f„), 

through the variable transformations, 

q n =\og^-, (A-3) 

Jn 
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In eqs. (A-2), g and g are the auxiliary variables and D is the bilinear differential operator defined by 

D k x f-g=(d x -d y ) k f(x)g(y)\ y=x . 

Letting 

fn = x„(n), gn = x„-i(n), gn = x„+i(n), (A-6) 

we have three bilinear equations, the TL equation, the Backlund transformation (BT) for TL and the 
DTL equation, 

D% (*) ■ T„ (*) = 2(T n+1 {k)x n . i (*) - X n (*)l n (*) ) , 

(oD* - 1 • x„_ i (* - 1 ) + x n (k - 1 )T„_ i (*) = 0, (A-7) 
x n (k+\)x n (k-\)-x n (k)x n (k) 

= a 2 (x„+i (k + 1 )T„_ i (k - 1 ) - X„ (fc)T n (*) ) . 



Appendix B: Proof of Eqs. (5.11)-(5.13) 



Let us first consider the Casorati determinant D = det (pj i A i + qj B;). Since each row is 

l<i,j<N ' ' 

sum of two vectors, the determinant is given by the sum of 2 N terms, 



D = 



Bi qi B! ••• q N l ~ x B l 

B 2 q2^2 ••• q% l B 2 

B N qN^N ■■■ In B n 

B 1 •• 



N 

+ 1 



Bi qi Bi ■■■ 



Bn qpjBjy ■■■ q N 1 Bn 



+ I 

l<Il<!2<iV 



Aii Pi l A{ l ■ ■■ pf { 1 Aj l 



Ai 2 Pi 2 Aj 2 



Pi 2 A n 



Bn QnBn ■■■ 1 Bn 

N N 



+ ■■■ + 



M P\M 
A 2 p 2 A 2 

An PnAn 



P^An 



A(9i,92,--- ,m)Y[ Bi+ L A (<7i>--- ,Ph,--- ,qN) A h Y[ Bi 



1=1 1 1 = 1 



N 



+ L A (iu--- ,Ph,--- ,Pi 2 r-' ,qN)A h A i2 Y[ B t + - 

\<h<h<N 



i=l 
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N 



£ A(q u --- ,p h ,--- ,p h ,--- ,p iM ,--- ,q N )A h Ai 2 ---A iM Y[ B i 

l<h<i2<~<iM<N 



i=l 



N 



H \-A(pi,p 2 ,--- ,PN)Y[ A i 



i=\ 



N 



N 



= A(qi,q 2 ,--- ,qN)Y[ B i 1+ L 



A(?i,-" ,Ph,--- ,lN)A h 



i=\ 



+ y ,Pi 2 ,-" ^jv) A h A ii { 

i<h<i 2 <N A (qu--- ,1n) B h B h 



+ I 

1<!1<!2<-<'V<A' 



A(gi,-" ,Ph,-" ,Pi 2 r" ,PiM>"- ,4n) yrAiy 



where A is the Vandermonde determinant defined by 



A(x u x 2 ,--- ,x N ) = Yl ( x i~ x j)- 

N>i>j>\ 



By rewriting as 



N 



Pi - qk 



1 if « - ik 

each term of the summation is given by 

A(?i,-" ,Ph,--- ,*v) .jrdf» 

V=l 



A(?i,- 



A(q h ,q h ,--- ,q iM 

MPhiPhi--- ,Pi U 



A([q h ,q h ,--- ,q iM 
thus we obtain 



M AT 

n n 

i¥>hfa>-m 
M M 



Piy-Vi 



M I N 

n n 



<7i'v - % I 



V=l j=l <7<v <7/ J v=l I *=1 Ph Qk J B'i 

His, 



nn 



9iv - 9«, 



n 



(PL-PiJiWu-liJ 



M 



\l=l»=lPh lip } v=l ^iy \1<//<V<M 0% 1iv)(lip Ph) J v=l 



V ;=1 / M=0 \<h<i 2 <-<i M <N 



( i-r K-Av)(^-^v) \ AA, V 
11 f„. _ V„. _ „. ~\ I 1 1 R' ' 



\l</i<V<M (/% - #v ) - Pis ) / v=l B l 

The prefactor of above summation gives the gauge factor in eq. (5.10), and from the summation part, 
F n , G n and H n in eqs. (5. 1 1)-(5. 13) are derived by taking 

qi = ^jr, Pi=a + \J a 2 -\r u |r ; | = 1, 
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A ; = />?exp&), 3 = 4nexp(-Si). fbrF„, 

Pi 

P n+1 1 

A,- = expfe), = exp(-K), for G„, 

A i = tf- 1 (l-«p i )exp&), s;- = ^^exp(-^), for/*,. 

Pi 
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